Coulomb Charging at Large Conduction 



Xiaohui Wang and Hermann Grabert 
Fakultdt fiir Physik, Universitdt Freiburg, Hermann- Herder-Strajie 3, 79104 Freiburg, Germany 

(February 1, 2008) 

We discuss the suppression of Coulomb charging effects on a small metallic island coupled to 
an electrode by a tunnel junction. At high temperatures the quantum corrections to the classical 
charging energy Ec = e^/2C, where C is the island capacitance, are evaluated. At low temperatures 
the large quantum fluctuations of the island charge cause a strong reduction of the effective Ec 
which is determined explicitly in the limit of a large tunneling conductance. 



, In recent years Caulomb blockade effects in metallic nanostructures have been studied extensively both theoretically 
' and experimentallytJ. These effects arise in systems with small metallic islands separated by tunnel barriers. For the 
T— I ', case of weak tunneling addressed mostly so farl3 the transfer of charges can be described in terms of a tunneling 
^ ■ Hamiltonian treated perturbatively. On the other hand, for strong tunneling one expects a washout of charging 
^ [ effects by large quantum fluctuations of the charge. For a quantum dot coupled by a point coijtact to a lead electrode, 
. the dependence of charging effects on the transmission coefficient was studied in recent workH. However, in this case 
. ' there are only few transport channels. 
^ . The behavior of metallic tunnel junctions with many channels can conveniently be described by a path integral over 
I ' fluctuations of the phase conjugate tostie island chargea. Based on this approach the suppression of Coulomb effects 
' was studied earlier by several authorga~l3. Most of these papers have considered single junctions while it become clear 
(Sj laterO that in this case charging effects are usually suppressed by the coupling to electromagnetic modes. On the 
Cn| other hand, charge quantization-oiL-sinall metallic islands is predominantly affected by fluctuations of charge carriers 

■ across adjacent tunnel junctionsEj tj. 

' The simplest device displaying these effects is the so-called single electron box (SEB) which is formed by a metallic 
islandpbetween by a tunnel junction and a gate capacitor (cf. Fig. 1). The external electrodes are biased by a voltage 
sourcalj. Recently, Falci et al.u have examined the low temperature behavior of the SEB by means of renormalization 
group techniques. While powerful in predicting the main parameter dependence of the effective charging energy in the 
' low temperature range, this method does not allow for a determination of the precise relation between low temperature 
Ch ' parameters and quantities measurable in the classical regime. On the other hand, experiments cover both the range 
I 1 of classical and quantum mechanical charge fluctuations and one would like to examine the low temperature behavior 
■ O ' of devices in terms of parameters measured in the classical regime. 

C , In this article we present a new method to evaluate the path integral in phase representation. At low temperatures 
the path integral is dominated by "sluggon" trajectories describing phase creep by 2tt across the junction. The 
contribution of these trajectories can be evaluated explicitly in the limit of low temperatures and large tunneling 
conductance. As a main result we find that the dominant effect of charge fluctuations can be absorbed in an effective 
charging energy. 

■ We start from the partition function of the SEB at inverse temperature (3 = l/ksT which may be written as 
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where the integral is over all paths of the phase ip{t) in the interval —(3/2 < r < /3/2 with if {(3/2) = if {—(3/2) modulus 
27r. The action S'box[¥'] — Sc[f \ + St[f \ contains two parts describing charging of the island and tunneling across the 
junction. The effect of the Coulomb energy is contained in 
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where Ec = e^/2C is the single electron charging energy, in which the island capacitance C is the sum of the 
capacitance Ct of the tunnel junction and the gate capacitance Cg. ricx — CgU/e is a dimensionless measure of the 
voltage U applied via the gate capacitance (cf. Fig. 1). The second part of the actiorD 
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describes tunneling. The Fourier transform of a(r) reads a; = — at|w/|/47r for |a;;| <C D where the uji = 27rZ//3 are 
Matsubara frequencies, D is the electronic bandwidth, and at = Rx/Rt is the dimensionless conductance of the 
junction, where Rt is the tunneling resistance and Rk = h/e'^ the von-Klitzing resistance. 

Since the action is a periodic function of ip with period 2tt, the partition function may be written as a sum over 
"winding numbers" , 
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Here 



Zk = 
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where the action S[(p] ~ S'oM + St[(p], with So[(p] = y^p^^dTip^ / AE^. 

Due to the nonlinear, nonlocal interactions described by St^-p]^ the partition function cannot be evaluated exactly. 
Here we shall employ semiclassical methods. The classical paths (from 5S[^p\ ~ 0) satisfy 
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ap^the boundary condition for winding number k reads (p(/3/2) = Lp{—j3/2) + 27rfc. A trivial solution of this equation 



with the classical action 
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For an arbitrary path v?(t) = 'fiiW''') + (^(t) of winding number k fluctuating about the classical solution, the second 
order variational action reads 
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dT'a{T - t') cosK(r - r')][6l(r) - 0{t')] 
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In terms of the Fourier transform 6i = 6\ + i6'", the action takes the form 

oo 

where for I <^ (3D 
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For very large / the eigenvalues approach ujf /2Ec independent of k. Upon normalization of the path integral the large 
I contributions cancel. It is now straightforward to show that 
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where Zq is the contribution of paths with winding number which is independent of ricx and 
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where u = at(3Ec/2TT'^ and k± = k + ^ zL ^y/Auk + u^. At high temperature this result is of the form of the classical 
result 

oo 

exp[-/3£;,(n - nex)'] (4) 
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yet with a renormalized charging energy 

E: = E,{1 - atl3Ej2Ti'' + 0[{l3E,f]}. 
The semiclassical evaluation of the partition function is well-behaved for all at if PEc ^ tt^. However, at low 

(k) 

temperatures the "time scale" j3 is large and the eigenvalues AJ for winding number k are nearly zero for small 
I < \k\. Thus, a simple semiclassical approximation with Gaussian fluctuations around the classical paths becomes 
obsolete. On the other hand, a vanishing eigenvalue indicates free motion of the minimal action trajectory in the 
direction of the corresponding eigenfunction. Hence, we expect that at low temperatures our system has a family of 
trajectories of almost the same action. 

Let us focus on the case on k = 1 where Lp{(3/2) ~ if{—(3/2) + 2tt and consider trajectories which gain a phase 
change of 27r within the imaginary time interval s < (3. For given s there is one trajectory with smallest action. 
This action is found to be almost independent of s as long as s ^ l/Ec- For s — (3 one recovers the classical path 

'P^^ ■ Hence, there is a family of trajectories that have nearly the same action. In contrast to the familiar instanton 
trajectories determining the low temperature behavior of multi-well systems, the trajectories considered here do not 
have a well-defined width, rather their action becomes smaller as the width grows. In view of their preference of a 
sluggish phase change, we shall refer to these trajectories as "sluggons" . 

In general, the sluggon trajectories cannot be found analytically. Here we consider the case at S> 1, l3Ec 3> 1 which 
allows for an analytic treatment. In this limit multi-sluggon trajectories with many phase changes of 27r dominate 
the path integral. The imaginary time interval /? may then be decomposed into segments of length s each containing 
a single phase change of 27r. Hence, a typical sluggon is constant outside a time interval s <C /? and moves within the 
interval s along the path of least action. Now, in the limit considsced, the kinetic part of the action may be treated 
perturbatively, and one finds from the equation of motion (|l|) thattj 

if>{T) = 2 arctan(fiT) + tt + 0(l/at, l/r2s), (5) 

is almost a minimal action path for arbitrary within the interval 1/s ^ <C -Ec- The sluggon action reads 

which is very weakly dependent on f2 as long as U, <^ Ec- 

The eigenvalue problem for fluctuations 9{t) about the sluggon trajectory (p{t) reads 

^"^^^ 2 f'^ dT'a{T ~ r') cos[^(r) - ^(r')][0„(r) - 0„(r')] + A„0„(r) = 0, 
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with the boundary condition 0„(— s/2) = 6'„(s/2). The eigenvalues A„ [n — 1,2,3,...) are discrete and doubly 
degenerate. In the limit a* » 1, s » l/fi » 1/-Ec, the low-frequency eigenvalues may be determined explicitly. We 
find 
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where Vn = 27rn/s. For n — 1^ the eigenvalue is nearly zero, and there are two linearly independent zero-modes that 
may be chosen as 
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and 
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These modes are associated with a variation of the shiggon width and the sluggon center, respectively. In the time 
interval — s/2 < r < s/2, the eigenvalues of the fluctuation modes about the trivial trajectory <^(t) = are simply 

^0 ^ . (^tVn 



2Ec 2tt 

The van Vleck determinant for fluctuations about the sluggon with the zero-modes omitted and normalized by the 
fluctuation determinant about (fir) — may be evaluated according to 
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which gives for atsEc ^ 1 
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The contribution to the partition function of paths with winding number k is the sum over all trajectories containing 
m -\- k sluggons and m anti-sluggons. Neglecting inter-sluggon interactions, one readily finds 

^ {2m + k)\ '^'^ n^+' , n^-^- , 2Xe-"'/2 

The integrations over the amplitudes of the two zero-modes of each sluggon have been replaced by integrations 
over the two "collective coordinates" T2p-i and T2p of the sluggon, where (T2p-i + T2p)/2 is the sluggon center and 
T2p — T2p~i = ■n/VLp the sluggon width. The factor 2/(T2p — T2p-i) is the Jacobian resulting from this transformation 
and Tc is the smallest sluggon size which is of order 1/ Ec- 

To estimate the effect of inter-sluggon interactions, we first note that the two sluggon contribution to Z2 reads 

^ = 2(/3X)2e-"* ln^(/3/2Tc). 

This result will be modified by inter-sluggon interactions. For two sluggons with centers at ti and T2, and widths 
7r/r2i and tt/^2, respectively, the asymptotic form of the interaction for large separation is —at/^i0.2{Ti — T2)'^ ■ When 

this is taken into account, we obtain a correction to Z^"* proportional to fp' ^ hence smaller than by a factor of 
order ln~^ (/3/rc). On the other hand, for small separation the action of two sluggons including interaction is bounded 
by the action of a 47r-sluggon, which is of order /31n(/3/Tc). Hence, both for small and large separation one finds a 
negligible correction if /3 is sufficiently large. Likewise, n non-interacting sluggons give a leading order contribution 
to Z proportional to (/31n/3)" while interactions give only corrections of smaller order for large /3. This indicates that 
the approximations made in deriving are indeed reasonable at low temperature. 
Now, from (nh the contribution for winding number k reads 
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and the partition function may be written in the form (|^) with Ck = Zi^/Zq. This result may be evaluated further 
for very low temperatures. Let us rewrite (|^) as 

00 

Zfc=^e/-W, (8) 

where 



fk{x) = X In [2/3Kexp{-at/2) ln{(3/Tcx)] - In [r(2;/2 + k/2 + l)T{x/2 ~ k/2 + 1)] . 

In the limit at 3> 1, PE,, 2> 1, this function has a maximum at a;,„ = 2/3i^at exp(— Q;t/2)[1 — 6lTi{at)/at + 0{l/at)] 
provided k <C x^m and the distribution exp[/j(a;)] becomes Gaussian with a width of order (2x^)^/2^ Thus the 
relevant contributions to result from sluggon configurations with sluggon numbers within Xm ~ (2a;m)^^^ and 
Xm + (2xm)^^'^ . The mean distance between sluggons is 13/xm, which is of order a^^ e-xp{at /2) / Ec- This is much 
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larger than Tc for large at, but much smaller than /3 at sufficiently low temperatures. Hence, the conditions under 
which sluggon trajectories of the form (||) dominate are satisfied. |— ■ 

The sum in may now be evaluated by means of the Euler-Maclaurin formulall3 yielding 

Ck = exp {-7T^ky/3E: + 0[ky{l3E:f]) , (9) 

where 

E*JE, = 2a3e-"'/2[l + 0(ln(at)M)]. (10) 

For riex not near 1/2 modulus 1, terms with large winding number are irrelevant and the approximation (^) suffices 
to determine the partition function, which is found to be of the classical form (^) yet with the charging energy E* 
defined in (pO|). Hence the quantum corrections may be absorbed in an effective charging energy. The exponential 
dependence of E* on at is in accordance with the renormalization group analysisH. The approach presented here 
determineS|-the preexponential factor explicitly for large tunneling conductance. The result (|lO| ) differs from earlier 
predictionsti for the effective charging energy of single junctions. 

In summary we have exploited semiclassical methods to evaluate the partition function of the SEB. We have obtained 
the leading order quantum corrections to the effective charging energy for arbitrary tunneling conductance at at high 
temperatures. At low temperatures the renormalized charging energy was determin|ed in the limit of a large tunnel 
conductance. The limit of small tunneling conductance has been treated previouslyllS. Combining these results one 
sees that there remains a lacuna at intermediate at where further work is needed. The result ([l0| ) can be extended to 
smaller values of at by avoiding the Euler-Maclaurin expansion and evaluating numerically, however, this is not 
sufficient to close the gap between the available small and large conductance results. 
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